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A  LIAPUNOV  FUNCTIONAL  FOR 
LINEAR  VOLTERRA  INTEGRODIFFERENTIAL  EQUATIONS 

D. L.Abrahamson  and  E.F.  Infante 


ABSTRACT 

Liapunov  functionals  of  quadratic  form  have  been  used  extensively 
for  the  study  of  the  stability  properties  of  linear  ordinary,  functional 
and  partial  differential  equations.  In  this  paper,  a  quadratic  functional  V 
is  constructed  for  the  linear  Volterra  integrodifferential  equation  . 

x(t)  =  Ax(t)  +  [  B(t-T)x(t)dT  ,  t  >  t  , 

J0 

x(t)  =  f(t) ,  0  <  t  £  tQ  . 

This  functional,  and  its  derivative  is  more  general  than  previously 

i 

constructed  ones  and  still  retains  desirable  computational  qualities;  moreover 
it  represents  a  natural  generalization  of  the  Liapunov  function  for  ordinary 
differential  equations.  The  method  of  construction  used  suggests  functionals 
which  are  useful  for  more  general  equations. 
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1 .  Introduction 


For  ordinary  differential  equations  of  the  form  x(t)=Ax(t),  x(t)  £  IR 

it  is  well  known  that  stability,  perturbation  and  asymptotic  properties  can 

T 

be  obtained  through  the  use  of  Liapunov  functionals  of  the  form  V=x  Mx. 

Indeed,  if  A  is  asymptotically  stable  a  positive  definite  symmetric  matrix 

M  can  be  obtained  as  the  solution  of  the  Liapunov  algebraic  equation 
T 

A  M  +  MA  =  -W  for  any  positive  definite  matrix  W;  moreover,  one  has  the 


familiar  result  that  M  = 


T 

Ac  Ac  •  T 

eA  we  ds  and  that  V  =  -xwx.  Since  V 


and  V  are  quadratic  forms,  it  is  easy  to  obtain  an  estimate  of  the  type 
V(t)  £  -26v(t)  for  some  positive  number  <$  and  therefore  to  obtain  a 
decay  estimate  via  the  Gronwall  inequality.  It  is  well  known  that  this 
estimate  can  be  made  arbitrarily  close  to  the  best  possible  estimate;  further- 

i  /  /» 

more,  V  '  is  a  norm  equivalent  to  the  Euclidian  norm  in  IR.  The  simplicity, 
computability  and  sharpness  of  these  results  are  central  to  a  number  of  useful 
techniques  in  asymptotic  and  perturbation  analysis  of  ordinary  differential 
equations . 

In  this  paper  we  present  a  quadratic  Liapunov  functional  for  a  linear 
Volterra  integrodifferential  equation.  Because  the  function  space  involved 
is  not  IRn,  it  is  not  clear  what  quadratic  functionals  should  be  chosen  that 
will  lead  to  sharp  results  yet  remain  easily  computable.  We  follow  the  methodology 
suggested  by  the  ordinary  differential  equation  case  and  choose  what  seems  an 
appropriate  quadratic  form  V  for  the  functional.  Computation  of  the  rate  of 
change  of  the  functional  along  the  solutions  of  the  equation  yields  another 
functional  V;  for  the  sake  of  computability  and  comparability,  conditions  are 
imposed  on  V,  which  are  then  reflected  in  the  form  V.  In  the  course  of  this 


.>■>  r  iSy  i-Afrv  i 


analysis  an  auxiliary  equation,  the  natural  generalization  of  the  Liapunov 
algebraic  equation,  arises  and  is  analyzed.  The  end  result  is  a  theorem  that 
states  the  existence  of  Liapunov  functionals  of  the  desired  form  which  give  sharp 
estimates  of  decay  rates. 

We  consider  the  equation 

rt 

(1.1)  x(t)  =  A(t)x(t)  +  K(t,T)x(T)dT,  t  >  t  >  0,  x(t) €  Fn 

Jq  U 

with  x(t)  =  f(t)  for  f  €  C[0,tp].  Miller  [1]  and  others  have  shown  the 
desirability  and  usefulness  of  considering  such  equations  with  t^  >  0.  Denoting 
by  x(*,tg,f)  the  solution  of  this  equation  with  initial  data  f  €  CfO.t^],  we 
say  thnt  the  equation  is  stable  if  for  every  tQ  0  and  e  >  0  there  exists 

a  <5  =  S(t0,£)  >  0  such  that  |f(t)|  <  0  £  t  £  t  ,  implies  jx(t,t  ,f)|  <  £ 

for  t  _>  t^;  the  equation  is  said  to  be  uniformly  stable  if  it  is  stable  and 
6  =  6(e)  can  be  chosen  independent  of  t  _>  0. 

For  the  simpler  equation 

.  n 

(1.2)  x (t)  =  Ax (t)  +  B(t-t)x(t)dT,  t  >_  t  >  0,  x (t )  £  K  , 

J0  V 

where  A  is  a  constant  matrix  and  the  convolution  kernel  B  is  integrable, 
characterizations  of  stability  properties  have  been  obtained  by  Miller  [1]  and 
Grossman  and  Miller  [2].  Seifert  [3,4]  and  Grimmer  and  Seifert  [5]  have  studied 

this  equation  through  the  use  of  the  very  simple  quadratic  Liapunov  functional 

T  T 

x  Jx,  under  the  assumption  that  A  is  asymptotically  stable  and  that  A  J+JA 

is  negative  definite.  Burton  [6,7]  has  studied  equations  (1.1)  and  (1.2)  using 

more  general  functionals,  generally  not  quadratic  ones,  functionals  which  include 

the  history  of  the  solution.  The  construction  of  quadratic  Liapunov  functionals 

for  linear  problems  in  Hilbert  spaces  has  been  pursued  extensively  by  Infante 


-s- 


and  others  [8,9,10,11].  This  latter  approach  has  motivated  this  investigation; 
indeed  the  results  here  generalize  those  in  [9].  In  Section  2  we  construct  a 
quadractic  Liapunov  functional  for  equation  (1.2);  in  Section  3  we  exhibit  a 
functional  which  is  useful  for  the  study  of  equation  (1.1)  and  present  an  example 
of  its  use. 

We  wish  to  thank  our  colleague  Professor  John  Mallet-Paret  for  some  very 
useful  suggestions  which  have  been  incorporated  in  this  paper. 
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2.  Equations  with  Kernels  of  Convolution  Type 


Consider  the  equation 

(2.1)  x(t)  =  Ax(t)  +  (  B(t-T)x(i)dT,  t>  t  >  0,  x(t)£Rn, 

J 0  u 

where  A  and  B  are  n  x  n  matrices  with  A  constant  and  B  locally  integrable. 
For  each  t  >  tQ  we  consider  the  pair  (xt(0),xt)  €  1R:1  x  Lo[-t,0],  where 
xt(0)  =  x(t+0)  for  -t  <  0  <  0.  Associated  with  this  equation  is  the  equation 

(2.2)  y(t)  =  (A+  6l)y (t)  +  j  el5(-t‘T-)B(t-T)y(T;di,  t  >_  tQ  >_  0, 

where  6  € R ,  which  is  obtained  from  (2.1)  through  the  change  of  variables 
y(t)  =  e6tx(t). 

In  the  sequel,  we  shall  use  a  representation  of  the  solutions  of  (2.1) 
similar  to  those  given  by  Miller  [12]  and  Hale  [13].  Let  S  be  the  solution 
of  the  matrix  initial  value  problem 


(2.3) 


S(t)  =  S(t)A  + 
S(0)  =  I, 


ft 

0 


S(t-T)B(T)dT 


t  >  0 


and  let  S(0)=O  for  0  <  0.  Then  we  have 


Lemma  2.1.  Suppose  x(t)  is  a  solution  of  (2.1)  for  t  >  tQ.  Then  if 
t  >_  t  and  p  >  0 


(2.4) 


f  o 


x  (0)  =  S(p)x  (0)  +  S(p-v-0)B(v)x  (0)dvd0  . 

x  p  LtJ.Q 


The  matrix  function  S  (t)  =  eQ  S(t)  is,  of  course,  related  to  equation  (2.2). 

o 

The  quadratic  forms  to  be  used  as  candidates  for  Liapunov  functionals 


for  equation  (2.1)  are  of  the  type 
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(2.5) 


V(t,xt(0),xt)  =  x‘(0)Lxt(0)  +  j  x‘(8)H(e)xt(0)d0 


♦  2x^(0) 


N(0)xt(0)d0  + 


0  , 0 


t'-t 


x‘(0)P(0,n)xt(T^)dnd0 


where  L  and  H(0),  -t  <  0  <  0,  are  symmetric  positive  definite  matrices 
T 

and  P(0,n)  =  P  (rq0)  for  -t  <  0,n  £  0  .  This  quadratic  form  is  natural  for 
our  equation.  The  rate  of  change  of  V  along  the  solutions  of  (2.1)  for  t  >  t 
can  be  obtained,  after  a  lengthy  but  straightforward  computation,  as 


V(t,xt(0),xt)  =  x^(0) [ATL  +  LA  +  H(0)  +  N(0)  +  NT(0)]  xt(0) 


xJ(0)H'(0)xt (0)d0 


(2.6) 


+  2x^(0) 
0  f0 


[lb (- q)  +  A  N(q)  -  N’(n)  +  p (o, n) ]xt (n)dn 


-t 


-t  J-t 


x[(0)[BT(-0)N(n)  ♦  NT(0)B(-n)  -  -P|pq)-  -  ]xt(n)drxl0 


This  quadratic  form  is  of  the  same  type  as  (2.5).  Unfortunately,  at  this  level  of 
generality  it  is  difficult,  if  not  impossible,  to  obtain  simple  computable  inequalities 
between  V  and  V. 


In  [3, 4, 5, 6, 7]  this  difficulty  has  led  one  to  choose  all  or  several  of  the 

matrix  functions  H,  N  and  P  to  be  identically  zero,  yielding  rather  restrictive 

results.  The  analogy  to  the  case  of  the  simple  ordinary  differential  equation 
•  q- 

x  =  Ax  with  Liapunov  functional  V  =  x  Jx  is  to  pick  a  particularly  simple  form 

T  T 

J  =  J  ,  say  diagonal,  and  to  insist  that  the  form  E  =  - (A  J  +  JA)  be  positive 

definite.  It  is  well  known  that,  in  general,  the  stability  results  thus  obtained 

for  the  solutions  of  the  differential  equation  will  be  most  restrictive.  Our 

viewpoint  is  the  opposite  one.  We  choose  a  particularly  simple  form  for  the 

matrix  E  and  thus  obtain  the  matrix  J.  This  is  precisely  the  viewpoint 
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successfully  employed  in  [9],  and  which  we  follow  here.  More  specifically,  we 
endeavor  to  obtain  a  relationship  between  V  and  V  of  the  simple  form 
V  £  -kV. 

In  the  ordinary  differential  equation  case,  our  viewpoint  leads  to  an 

T 

analysis  of  the  matrix  equation  A  J  +  JA  =  -E.  It  is  well  known  that  if  A 
is  asymptotically  stable,  for  any  symmetric  E  this  equation  has  a  unique 
symmetric  solution  J  and  that  the  mapping  E  J  maps  positive  definite 
symmetric  matrices  into  (but  not  onto)  positive  definite  symmetric  matrices. 
The  following  analysis  is  a  generalization  of  this  well  known  result  to  the 
problem  at  hand. 


Lemma  2.2.  The  matrix  problem 


(2.7) 


Q'  (a)  =  CQ (ct)  + 
Q(a)  =  QT(-a), 


r 

J  n 


D(3)Q(a-|3)dB  ,  a  >  o 


a  £  li 


CQ(0)  +  Q(0)CT  +  f  [D(f3)QT(£)  +  Q((3)DT((3)  ]dg  =  X, 
J0 


V  • 

where  C  is  a  constant  matrix,  X  a  symmetric  constant  matrix  and  D(*)e  € 

LjCO.oo)  fl  L2(0,co)  for  some  Y>0,  has  no  more  than  one  exponentially  decaying 
solution. 


Proof :  It  suffices  to  prove  that  if  X=0  (2.7)  has  no  nontrivial  exponentially 

decaying  solutions.  Suppose  that  Q  is  a  solution  with  |Q(°0|  £  ke~V^ 


for  some  k,v  >  0.  Define  Q(s)  = 


pu  A 

e  S  Q(cx)da  .  Then  Q  is  analytic  in  s 


in  the  strip  |  Re  s  f  <  v.  Define  H(ct)  for  a  e  ]R  by 

(2.8) 

it  follows  that  H(a)  =  0  for  >  0  and  that 


H(tx)  =  Q'(a)  -  CQ(a)  -  |  D(B)Q(a-B)db  ; 

0 
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(2.9) 


H(a)  =  -CQ (a)  -  Q(u)CT  -  f  D(£)Q(^-£)d£ 

i  0 

.  oo 

-  Q(a+j3)DT(e)d£  ,  a  <  0. 

JQ 


Extend  the  matrix  function  D  by  D(8)  =  0  for  £  <  0  and  define  H(s) 

and  D(s)  analogously  to  the  definition  of  Q(s).  It  is  easily  seen  that  H 
is  analytic  and  bounded  in  an  open  half-plane  containing  {s|Re  s  <  0}  and  that 
D  is  analytic  and  bounded  in  an  open  half-plane  containing  (s|Re  s  j>  0}  .  From 
(2.8)  it  follows  that  there  is  a  strip  about  the  imaginary  axis  in  which 


(2.9)  H(s)  =  A (s)Q(s) ,  A(s)  =  si  -  C  -  D(s)  , 

and  since  Q(s)  =  QT(-s)  it  follows  that  H(s)A^(-s)  =  A(s)H^(-s)  in  this  strip. 

Consider  the  function  F(s)  =  H(s)AT(-s),  which  is  analytic  in  an  open  half  plane 

T 

containing  {s  |Res  <  0}  ;  the  relation  F(s)  =  F  (-s)  shows  tl.  t  F  is  an 
entire  function.  Moreover,  it  is  easily  seen  that  there  exist  positive  constants 
c^  and  c2  such  that  [F(s)|  <  c^  +  c _, ]  s  [  for  all  s  £  (L.  By  a  familiar  exten¬ 
sion  of  Liouville's  Theorem  it  follows  that  F(s)  =  F^  +  sF2  for  some  complex 
matrices  F^  and  F2.  Hence,  H(s)  =  (Fj  +  sF0)  (-si  -  CT  -  DT(-s))  1  for 
Re  s  <_  0.  From  (2.7)  and  (2.9)  it  follows  that  H' (a)  exists  and  decays  exponen- 

A 

tially.  Since  lim  s  H(s)  =  H(0)  and  since  X=0  implies  H(0)=0,  it  follows 
Re  s  -*■  - 00  /v 

that  Fj=F2=0.  Thus  H(s)  vanishes  identically,  which  implies  that  Q(s)  and 
Q(cx)  must  be  zero,  proving  the  lemma. 

Let  us  now  return  to  the  construction  of  the  Liapunov  functional  for 
equation  (2.1).  For  this  purpose,  assume  that  for  the  given  A,  B  and  <5 
equation  (2.2)  is  uniformly  stable  and  that  the  kernel  B(*)e^  6  L^  CO ,°°)  0  L2(0,°°). 

Let  e  >  0,  and  consider  the  problem,  obviously  related  to  equation  (2.7),  of  finding 


a  function  Q  within  the  class  of  exponentially  decaying  matrices  defined  by 

•  °0 

Q’  (oi)  =  [AT  +  (6-e)  I]Q(a)  +  e(6't)  CBT(8)Q(a-$)de,  a  >  o 

(2.10)  Q(a)  =  QT(-a),  a  €  R 

GO 

ATQ(0)  +  Q(0)A  +  Q  e(  '£)6[BT(6)QT(B)+Q(6)B(e)]de  =  K 

for  some  fixed  but  arbitrary  positive  definite  matrix  W,  to  be  selected  in  the 
sequel.  Lemma  2.2  above  guarantees  that,  among  exponentially  decaying  matrices, 
there  is  at  most  one  Q  satisfying  this  equation.  With  £(t)  =  e^  ^SCt), 

where  S  is  the  solution  of  equation  (2.3),  it  is  easily  checked  that 

,<*> 

(2.11)  Q(a)  =  S^u  +ct)WS^_e(u)du,  <*€  JR, 

Jo 

yields  a  solution  of  (2.10),  hence  the  unique  exponentially  decaying  solution  of 
this  equation.  It  is  this  matrix  function  Q  that  allows  us  to  construct  a  simple 
Liapunov  functional  of  the  form  (2.5).  Indeed,  in  that  equation,  let  the  matrices 
L,  H(0) ,  N(0 )  and  P(0,n)  be  defined  by 

(2.12)  L  =  Q(0)  +  M, 

(2.13)  H(0)  =  R  e2!-6'e)e 

(2.14)  N(6)  =  eC6"£)(v+0)Q(v+0)B(v)dv, 

(2.15)  P(0,r|)  =  f  f  eC6'e)(ri+ 0+U+v)BT(v)Q(n-e  +u-v)B(u)dudv, 

J-0  J-n 

where  the  positive  semidefinite  matrices  M  and  R  will  be  chosen  later.  A 
lengthy  computation  yields  that,  with  these  matrices,  equation  (2.6)  reduces  to 

• 

V(t,xt(0) ,xt)  =  -2(6-e)V(t,xt(0) ,xt) 

+  xJ(0)[-W  +  R  +  (AT  +  (6-t) I)M  ♦  M(A+(6-e)I)]x  (0) 


(2.16) 
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Moreover,  it  is  possible  to  expose  the  structure  of  the  functional  V  under 
these  conditions.  Substitution  of  equation  (2.11)  into  equations  (2 . 12) - (2 . 15) 
shows,  after  some  appropriate  rearrangements,  that  the  functional  V  is  of  the 
form 


(2.17)  V(t,xt(0),xt)  =  xJ(0)Mxt(0)  ♦  J°  xJ>)Rxt(0)e 


2(6-e)6 


d6 


(S(p)xt(0)  + 


0 


x  (S(p)xt(0)  ♦ 


-t 


-0 


S(p-v-0)B(v)xt(c‘)dvde)TWe2^~e^P 


fO 


/: 


S(p-v-0)B(v)xt(0)dvd6)dp  , 


J-t  -'-0 

which  can  be  rewritten,  on  account  of  Lemma  2.1,  as 


(2.18)  V(t,xt(0),xt)  =  x[(0)Mxt(0)  +  J  xJ(0)Rxt(0)« 

♦  f  xJ(P)Wx  (P)e2(<S*E)PdP  . 
Jn  r  c 


2(6-£)0 


d0 


Three  points  seem  worthy  of  some  emphasis  at  this  point.  The  first  is  that 
whereas  equation  (2.5)  is  a  very  general  form,  it  is  the  particular  restriction 
of  this  form  given  by  (2. 11)- (2. 15)  that  yields  the  simple,  manageable  results 
expressed  by  equations  (2.16)  and  (2.18).  The  second  remark  is  that  the  simplest 
computable  result  for  (2.16)  is  obtained  for  the  case  6  -£  =  0,  M=R=0,  in  which 
case  the  sign  of  V  is  strictly  given  by  W;  but  in  this  case,  all  that  one 
can  say  is  that  V  is  positive  and  difficult  to  estimate.  Finally,  if  in  (2.18) 
one  chooses  W=0  ,  then  one  recovers  the  commonly  used  functionals  for  this 
problem,  but  then  note  that  (2.16)  is  difficult  if  not  impossible  to  estimate 
appropriately  unless  the  matrix  A  is  asymptotically  stable. 

The  emphasis  here  is  to  obtain  estimates  through  appropriate  use  of  the 
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matrices  W,  M  and  R.  First  note  that  with  W  >  0,  it  follows  from  (2.18)  that 
(2.19)  V(t,xt(0),xt)  >  x^(0)Mxt(0)  +  j  xJ(e)Rxt(0)e2(<S"£)0d0  . 

Let  now  KIM  =  (CTMC)^2  for  £  6  IR11  with  the  same  notation  for  the  corresponding 
matrix  norm.  Now,  let  e  >0  be  given.  Observe  that 


t  r°  f°  ? 

2x*  (0)  J  MB(-0)xt(0)de  <  2|xt(0)lM  (  j  |xt(3)|J 


e2(6-e)0de) 


1/2 


*  ( 


J-t 


r(-e) | 2  e'2(6'e)0de)1/2 


Bf>lM  e<6'e)’11  L2(0.«)(Clxt<°)lM  *  H.0tUt(e,l«  e2(S'E,ede) 


for  every  c  >  0.  Hence,  if  we  choose  c  =  e/||  |B(*)e^  C^||  and 

l|B(.)e(6-eHT2 


R  = 


2|x‘(6) 


M,  it  follows  from  the  form  of  (2.18)  that 
(0 


MB(-0)xt(0)d0|  <  2C  V(t,xt(0ixt)  . 


'-t 

Thus  with  this  choice  of  R,  and  M  an  arbitrary  positive  definite  matrix, 
(2.16)  yields 

V(t,xt(0),xt)  <  -2(6-  2e)V(t,xt(0)  ,xt) 

+  xJ(0)[-W  ♦  R  +  (AT+(6-e) I)M  +  M(A+(6-e)I)]xt(0) 

so  that  W  can  be  chosen  as  a  positive  definite  matrix  to  yield 


V(t,xt(0),xt)  <  -2(6-  2e)V(t,xt(0),xt). 


(2.20) 
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We  may  appeal  to  the  definitions  of  Q,  N,  and  P  to  show  that 
for  some  k  >  0  independent  of  t. 


(2.21)  V(t,xt(0),xt)  <  k^(0)Mxt(0)  +  |°  x[(6)Rxt(0)e 


2(6-e)0 


de 


Along  with(2.19\  this  establishes  just  in  what  sense  V  measures  the 

1/2 

size  of  the  solution.  For  each  f^xed  t,  V  (t,',‘)  is  equivalent 
to  the  norm  f  |  •  j  |  given  by 

(2.22)  ||(xt(0),xt)||2  =  |xt(0)|2  ♦  f°  |xt(9)|2e2(6'E)ed6  . 

•'-t 

1/2 

It  can  be  said  that  V  (t,-  , •)  is  equivalent  to  the  usual  norm 
on  J?n  x  I^f-t.O);  however,  the  constants  describing  this  relationship 

1/2 

behave  badly  as  t  grows,  whereas  the  constants  relating  V  (t,-,-)  to 
(2.22)  are  independent  of  t.  The  preceeding  results  are  summarized  in 
the  following 

•  ft 

Theorem  2.1.  Consider  the  equation  x(t)  =  Ax(t)  +  B(t-x)x(x)dx  and  the 

jo 

functional  V  given  by  (2.5).  Suppose  that  o  €  1R  is  such  that 

fi.  . 

B(*)e  €  L^(0,<»)  H  L^CO,00)  and  the  equation  y(t)  =  (A+6l)y(t) 

t 

♦  |  e^^tT^B(t-x)y(x)dx  is  uniformly  stable.  Then  for  any  e  >  0 
there  exist  continuous  matrices  N,  P,  and  H,  with  H  positive  definite, 
a  positive  definite  matrix  L,  and  *2  —  >  0  such  that 
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(2.23) 


< 


V«-£)0<bV/2<  V^2(t>;<t(0).xt) 


/ 

|°  K®)| 


vH-e)e<»y/2 


and  V<-2(6-2e)V. 

If  (5  >  0,  the  problem  is  exponentially  asymptotically  stable  and 
V  yields  the  estimate  |xt(0)j  <  c0e” for  a  constant  c^ 
depending  on  the  initial  data  on  [O.tg]  giving  rise  to  x. 


Remarks .  The  construction  of  a  Liapunov  functional  for  the  equation 

(2.24)  x(t)  =  Ax(t)  +  f  B(t-l)x(x)dt  , 

■'t-r 

where  0  <  r  <  <*>  and  B  e  L2(0,r),  is  carried  out  in  exactly  the  same 

manner  used  above.  The  result  is,  in  fact,  more  satisfactory.  We  regard 
(2.24)  as  an  autonomous  evolution  equation  in  the  Hilbert  space  JR  x  L2(-r,0). 

If  the  spectrum  of  the  generator  of  the  appropriate  semigroup  lies  in  the 
closed  half-plane  Re  s  <  <Sq,  then  for  any  eQ  >  0  we  can  construct  a 
quadratic  form  VQ  on  JRn  x  L2(-r,0)  such  that  Vg^2  is  equivalent 

to  the  usual  norm  and  Vq  <_  2^o+£0^0‘ 

The  approach  taken  here  may  be  of  more  general  use.  The  linearity 
and  appearance  of  the  equation  under  consideration  suggest  a  quadratic 
Liapunov  functional.  Above,  the  attempt  to  obtain  a  certain  appearance 
in  its  derivative  required  that  terms  in  (2.5)and  (2.6)be  comparable. 

In  general,  a  set  of  equations  will  arise  from  the  comparison  of 


V  to  V,  and 


various  tactics  are  available.  For  example,  dropping  the  requirement 


T 

P(9,n)  =  P  (p,9)  leads  to  a  functional  which  is  useful  for  an  equation 

whose  kernel  is  not  of  convolution  type.  This  remark  is  pursued  in  the  followin 

section. 

3.  More  General  Equations 

We  now  consider  the  equation  in  IRn 

(3.1)  x(t)  =  A(t)x(t)  +  [  K(t,T)x(x)dT, 

Jo 

where  A  and  K  are  continuous  n  x  n  matrices.  The  initial-value 
problem  to  be  considered  has  tQ  >  0  and  x  =  f  on  [0,t0],  f  £  C[0,tQ], 

with  (3.1) to  hold  in  t  tQ.  Again,  for  t  >  tQ  we  think  of  (xt(0),xt) 

as  lying  in  IR11*  L2(-t,0),  and  seek  a  quadratic  Liapunov  functional. 

Burton  [6,7]  and  others  [3,4,5]  have  studied  (3.1)by  picturing  it  as  a 
perturbed  version  of  an  ordinary  differential  equation,  and  we  do  so  as 
well.  The  theorem  presented  in  this  section  is  said  to  be  a  Razumikhin- 
type  result  because  its  proof  relies  only  on  the  description  of  certain 
solutions  at  the  times  those  solutions  leave  sets  which  are  of  interest. 

Before  presenting  the  theorem,  we  recall  a  simple  fact  from  linear  algebra 
which  can  be  found,  for  example,  in  [14].  If  D  and  E  are  real  symmetric 
n  *  n  matrics  with  E  positive  definite,  the  equation  det[D-\E]  =  0 
has  n  real  solutions  They  coincide  with  the  eigenvalues  of  DE_1 
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i 

t 


As  before,  for  E  >  0  and  5  £  !Rn  ,  we  write  |£|  =  (^ES)*^,  with  the 


same  notation  for  the  corresponding  matrix  norm  .  If  G  is  an  n  *  n 


matrix  with  real  eigenvalues,  let  Xmov(G)  denote  its  largest  eigenvalue. 

IQoa 


Theorem  3.1.  Suppose  M  is  a  positive  definite  n  *  n  matrix  and 


T  -1 

let  g(t)  =  ^ay(A  (t)  +  MA(t)M  ).  Suppose  that  there  exists  a  >  o 


such  that 


(3.2) 

for  all  t  >  0.  Then 


f0  DKcu>s)|«(exp0) 


dsdu  <  1  -  a 


(i)  (3.1) is  uniformly  stable  if  there  exists  N  <  such  that 
,T 


I  g  (  x)  d  x  £  N  whenever  0  £  t  £  T  <  +°°  ; 


(ii)  (3.  l)is  stable  if  lim  sup  f  g(T)dx  <  +»  . 

T-*”  JO 


Proof .  Let 


(3.3)  V  (t,x(t),x)  =  xT(t)Mx(t)  -  2[  f  xT(u)MK(u,s)exp|  [  g)  x(s)dsdu. 

JO  J  0 


Then  along  solutions  of  (3.1) 


(£•) 


(3.4)  V1(t,x(t),x)  =  xT(t)[AT(t)M  +  MA(t)  -  g(t)M]x(t)  +  gftJVjtt.xW.x), 


Thus  Vj(t,x(t) ,x)  <  g(t)V1(t,x(t) ,x) .  So,  for  t>  tQ, 

(3-5)  V1(t,x(t),x)  £  Vj^.fCt^.fJexpjJg^dx  , 


-17- 


where  f  is  the  initial  data  for  x.  If  (3.1)were  not  uniformly 

00 

stable,  we  could  find  sequences  as  follows:  {t^}^^  >  0, 

{Vk=l'  fk  6  C[0,tk^  ^UP  lfkCt^M-  lf  and  {Vk=l’  Tk  —  V  SUch 

that  k  and  |Xk(t)|M  <  k  for  t  <  T^,  where  x^  is  the 

solution  of  (27)  arising  from  initial  data  f^  on  [0,t^].  For  any 

k^l, 

Vl<VW’Xk>  i  I  WlJ 

rT,.  cu 


>  <xk2. 


Also, 


v!tvW'Viiwi2 

+  2  I C)k  lo  *fk(u)lM^U's)lM  eXp^|uk  fkCs)lMdsdu 
<  1  +  2  |  K(u,s)|m  expfj  k  gj  dsdu 


<  2  -  a  . 

From  (3.5)applied  to  these  inequalities, 
(3.6) 


ak2  <  (2-a  )exp  (  k  g(x)dx 


Since  (3.6)must  hold  for  every  k,  lim  g(x)dt  =  This  proves  (i) 

k~  lk 


i: 

ji 


Munuiiai 


p 
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To  prove  (ii)  we  follow  the  same  argument  except  that  all  the  t^ 
coincide  at  some  tQ  >  0.  The  same  chains  of  inequalities  lead  us  to 
this  version  of  (3.6): 


Otk^  <:  (2-a)exp  [  *Cg(t)dt. 

•>  t. 


This  implies  that  lim  sup 

t-*°° 


g(l)di  =  ■*»  and  thus  establishes  (ii) . 


Remarks .  The  formulas  in  (3. 3) and  (3.4)  resemble  those  describing  V 
and  V  in  Section  2.  In  fact,  is  motivated  by  the  construction 

there.  If  K(t,x)  =  (t-T)  and  A(t)  =  A.,  then  g(t)  = 

and  the  final  term  in  (3.3)  takes  the  form 


0  f0 


-tJ  -t 


xj  (a  )MB1  (0  -n)  e"g  l6  xt  (n)  dnde 


after  a  change  of  variables.  It  may  be  rewritten  as 
r°  r°  T 

j  j  xt(0)P1(0,n)xt(ri)dndW 

with  P1(0,n)  =  -MB^O-pie  gl®  for  -t  <  q  <  0  <_  0  and  P^n.Q)  =  Pj(^,n)- 


Then 


3px  (t3 ,  n)  aP^e.n) 
30  +  :  an 


=  -g1P1(e,n) 


so  that  Vj  fits  the  term-by-term  comparison  scheme  used  before;  the 
kernel  corresponding  to  N  of  Section  2  vanishes. 

The  simple  appearance  of  and  does  not  even  depend  on  the 


linearity  of  (3.1).  However,  we  have  not  as  yet  thoroughly  exploited  this. 


Functionals  similar  to  appear  in  Miller  [12]  for  the  case  of  a  scalar 

equation. 

The  functional  may  be  useful  in  analyzing  equations  whose  kernels 

do  not  decay  exponentially,  provided  the  equation  can  be  viewed  as  a 
perturbation  of  a  suitably  stable  ordinary  differential  equation.  The 
following  example  illustrates  this.  Consider  (3. 1)  where  n  =  2, 

| K(t,T)  |  £  k0min(l,Ct-T)-2),  and 

-1  3(t)  1 

0  -ij’ 

where  | 3(t) |  £  2-$0farsome  fixed  • >  0.  We  apply  Theorem  3  . 1  with 
M  =  I.  Then  g(t)  =  -2  +  | 3(t) |  £  -$Q.  Moreover, 


and  the  equation  is  uniformly  stable  if  4k^  £  (3g. 
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